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Abstract 

The results of a model for meson-meson scattering are studied. The model is shown to 
be capable of on the one hand reproducing the scattering data, while on the other hand a 
quark-antiquark confinement spectrum can be determined. 

It is concluded that adopting the model's formulation of the transition matrix elements 
for data analysis , it may serve as a link between experiment and quenched lattice calcula- 
tions. 



1 Introduction 

In the sixties, early seventies, mesons and baryons were studied in terms of models for confinement 
of the newly invented quarks |l|2j . The basic idea was that mesons and baryons could be described 
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by permanently confined quarks and/or antiquarks. In Fig. we give an artist's impression 
of the resulting mass spectrum for mesons. 
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Figure 1: Mass spectrum for a permanently closed system without (a), and including short 
distance effects perturbatively (b). The parameters are taken from Ref. [H], where harmonic- 
oscillator confinement is employed. 

Then, in the seventies, early eighties, short-range effects had to overcome the small disagree- 
ments between Nature and confinement models, as schematically depicted in Fig. (QJo). Later, 
hadronic decay was implemented, with effects on the meson masses as shown in Fig. |21 



2.0- 



1.5- 



1.0- 



0..-J- 



mass 






— 2D 




— 3S 






— IS 




— IS 


P-wave scattering 




perturbative 





(a) 



GeV mass 



2.0- 



1.0- 



0.5- 




S-w&ve scattering 



perturbative 



(b) 



Figure 2: Mass splittings for a permanently closed system when hadronic decay is perturbatively 
taken into account, for P wave (a), and for 5* wave (b). 

However, it was all the time overlooked how those mesons are actually produced in experiment. 
In particular, it was taken for granted that mesonic states below a particular meson-meson 
threshold do not feel the presence of that two-meson channel. Effects of coupling were assumed 
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to abruptly vanish or show up, depending on whether the state was just above or just below 
threshold. For states far below threshold, the whole idea of implementing two-meson channels in 
the confinement models was considered inopportune. Furthermore, little attention was paid to 
the fact that most mesons appear as resonances, and that mesonic resonance widths are usually 
not small. For example, a p meson has a lifetime compatible with its internal frequency, implying 
that it should certainly not be considered a more or less stable particle. Even nowadays this is 
not fully accepted or even understood. 

The procedure that mimics the experimental situation is called unitarization at present. This 
name suggests that all is accounted for, which is of course an illusion. In practice, one is more 
than satisfied if the "most important" effects are included in the resulting meson model. 

Also in the seventies one became aware of phenomena in 5-wave meson-meson scattering 
which could not be easily handled by any of the models. Hence, still in the spirit of permanent 
confinement, quark configurations were invented for mesons and baryons other [1J|I3II>] than the 
usual qq and q 3 . 

However, things really fall in place when full scattering properties are determined for meson- 
meson, or meson-baryon, processes, since then resonances come out automatically without any 
need to worry about their composition. 

We observe the following: 

• The higher radial excitations in cc and bb vector states are almost equally spaced in mass. 
Why the ground states for those systems come out much lower is well explained [Zj. 

• The level splittings for the S and D cc and bb vector states follow naturally [7j. 

• The phase shifts for Kir and irn in P wave |3], and in S wave jH], as well as the scattering 
lengths are reproduced. 

• The J p = + cs experimental results are explained [TU] . 

• The J p = 1 + cn (n for nonstrange, either up or down) and cs phenomena come out well [TT] . 

• It fully explains the light scalar meson nonet jH]. 

No extra forces and/or configurations are necessary to explain such an amount of very different 
data, with one set of parameters for all mesons. 

We may thus conclude that unitarization works well. Nevertheless, our model is far from 
perfect for many reasons that are easy to understand. First, for the transition potential, which 
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should have been derived from the theory of strong interactions, we make an as simple as pos- 
sible choice, and then it is still further simplified ^2]- Furthermore, also harmonic-oscillator 
confinement does not have much more justification than the equal-level spacings for cc and bb 
(and possibly also nn and ss) vector states. Hence, our choices for both, the confinement and 
unitarization ingredients, are taken as simple as possible. The fact that it works even under such 
extreme simplifications clearly indicates that, in a more elaborate approach, the unitarization 
scheme has a good chance to survive all tests. 

2 Unitarization 

In order to understand our approach to unitarization, one may imagine a huge scattering or tran- 
sition matrix, describing all possible scattering of meson pairs. Elastic processes appear on its 
diagonal, whereas off-diagonal matrix elements describe inelastic processes. Conservation laws 
predict that the vast majority of off-diagonal matrix elements vanish. The remaining nonvanish- 
ing matrix elements can be regrouped into smaller submatrices of meson-meson channels which 
under the given conservation laws are allowed to communicate. Hence, if we for example study 
the J p = + cs system, then we consider the part of the T matrix which describes the S"-wave 
scattering of all meson pairs that through OZI-allowed processes couple to cs. In practice, we 
limit ourselves to a few channels which we "believe" are relevant to the energy domain under 
study. In the following we shall concentrate on just one scattering channel. But the results are 
equally valid when more channels are taken into account. 

Let us now consider an arbitrary diagonal element of the above-discussed transition matrix in 
a specific partial wave, which we denote by T#. It describes the elastic scattering of two mesons, 
assumed not to couple to anything else but the confinement system having the corresponding 
quantum numbers. We assume here that confinement yields an infinite spectrum. The function 
Ti must be analytic in the total invariant mass E = yfs of the two-meson system. Moreover, it 
depends on the meson masses Mi and M 2 , the intensity A of the coupling between the two-meson 
system and the confined quark-antiquark system, and finally on the confinement spectrum E n 
(n = 0, 1, 2, . . .) and other details T q q of the qq system. In Ref. [T2"] . the interested reader may 
find a rather general and complete expression for 7}. Summarizing, we conclude 

T e = T t (y/Z, X, M 1 ,M 2 ;{E n \n = 0, 1, 2, ...},...) . (1) 

In Fig.|3]we have depicted a possible spectrum for confinement. The spectrum corresponds to 
A = 0, i.e., having discrete energy eigenvalues of the confinement operator under the assumption 
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that the qq system is decoupled from the meson-meson continuum. 
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Figure 3: The spectrum of confinement. 

When the confinement states are weakly coupled to the meson- meson continuum (A small), 
we observe narrow resonances in the two-meson scattering cross section a, about at the energy 
eigenvalues of confinement. This is depicted in Fig. 
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Figure 4: Confinement spectrum as observed in elastic meson-meson scattering, for the case that 
the confinement states are only weakly coupled to the meson-meson continuum. 

The cross section in Fig. H] is calculated by the use of the transition matrix element Ti of 
Eq. an explicit expression of which can be found in Ref. [12]. Now, since Tg is an analytic 
function in the total invariant two-meson mass E = y/s, one may study its singularity structure 
for complex E. This can be done numerically. In Fig. we show the poles associated with the 
resonances in the cross section depicted in Fig. 0] 
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Figure 5: The scattering-matrix poles associated with the resonances shown in Fig. 01 Threshold 
is at y/s = Mi + M 2 . The confinement energy eigenvalues are indicated by small dots on the 
real-energy axis. 
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The real parts of the poles are close to the energy eigenvalues of the confinement spectrum, 
while the imaginary parts are relatively small as the resonances are narrow. 

In the following we study how the scattering-matrix poles move through the complex-energy 
plane, when the coupling A between the meson-meson and confinement systems is increased, in 
particular near threshold. 

But first, let us discuss how a pole of the scattering matrix behaves when the confinement 
spectrum (A = 0) has a state which is below threshold. In our model, these poles shift to even 
lower energies when A is increased [7j , but remain on the real-energy axis. Such poles correspond 
to bound states of the coupled system. In this case, both the two-meson and the qq components 
of the wave function describe bound systems. The corresponding states are thus mixtures of 
bound two-meson and qq states. 

3 Behavior at threshold 

When the ground state of the confinement spectrum is near a non-^-wave threshold, then the 
associated T-matrix pole passes through threshold when A is increased ^T]. This phenomenon 
is depicted in Fig. |H1 



Figure 6: Scattering-matrix pole position of the ground-state resonance, as a function of the 
coupling constant A, for P- and higher-wave scattering. Some values of A are indicated along the 
curve. 

We observe that, while the real part of the pole is continuously decreasing when A is increased, 
the imaginary part first increases (in modulus), then decreases, vanishes right at threshold, and 
remains zero when A is still further increased. For small values of A, the displacements of the 
scattering-matrix poles in the complex-energy plane are almost perturbative, i.e., linear in A 2 . 
However, when A grows it becomes clear that the model incorporates highly nonperturbative 
effects. 

The pole displacement depicted in Fig |H1 represents two-meson scattering in P wave. For 
higher waves the behavior is similar. This can easily be understood by studying the effective- 
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range expansion near threshold. For P wave and higher waves, one finds that the imaginary part 
of the pole is quadratic in A when the real part of the linear momentum is small, which is near 
threshold. Consequently, at threshold it vanishes. 

However, S-wave poles behave very differently ^HIE]- They approach the real-energy axis 
perpendicularly well below threshold for increasing A. Then, for still larger values of A, they first 
move in the positive direction towards threshold, as virtual bound states. At threshold they turn 
back, now moving towards smaller energies as genuine bound states. This can more easily be 
understood from their movement in the complex- linear-momentum plane [ToTll3j . Experimentally, 
one can determine whether the pole represents a virtual or a real bound state from the sign of 
the scattering length. 

4 Poles from the quark-pair-creation cavity 

For small coupling, it is very clear how poles in the scattering amplitude are related to the energy 
eigenvalues of confinement, since their displacements are perturbative, hence small. However, for 
strong coupling pole displacements can be of the order of magnitude of the level splittings of the 
confinement spectrum. In thit case the situation is more complicated, and were it not for a model 
to trace the poles, the classification in terms of a confinement spectrum might turn impossible. 

But there is more: even when all poles are traced, starting from the small-A positions near the 
confinement spectrum, and ending up at the physical positions fitting the scattering data, there 
exist still more poles in the complex-energy plane. This means that the scattering amplitude 
which agrees with experiment contains more poles than just those stemming from the confinement 
spectrum. But let us first discuss where such extra poles happen to originate in the model. 

In our model, the two-meson system communicates with the confined qq system through 
OZI-allowed quark-pair creation. The shape of the corresponding potential has a maximum at 
relatively large distances (0.5-1.0 fm). This means that in the interior one has a small potential 
well which in principle could host bound states. For small coupling, the well is almost flat 
and thus of no consequence, giving rise to merely mathematical poles with very large negative 
imaginary parts, hence unobservable. However, for strong coupling these cavity poles will turn 
out to get mixed up with the poles originating in the confinement spectrum. It such a case, there 
are clearly observable effects in the scattering cross section. Under certain conditions, the poles 
even end up on the real-energy axis as bound states. In the past we referred to these poles as 
background poles, since for decreasing coupling they disappear into the background. Here, we 
shall stick to the term extra poles. 

The extra poles were first reported in Ref. [H|. Later, it was thought that they were the result 
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of pole doubling ^B] , a term which we adopted for a while. But now it has become clear from 
which mechanism they originate. 

For meson-meson scattering in S wave, it is more likely that certain structures in the cross 
section stem from the extra poles than for P and higher waves. The reason is that in .S-wave 
scattering the centrifugal barrier is absent. Consequently, the transition-potential well is deeper 
in the latter case. This is exactly what is observed in experiment. 

There are several parameters influencing the displacements of the extra poles, one of them be- 
ing the relative position of threshold with respect to the confinement spectrum. When threshold 
is far below the confinement-spectrum ground state, which is usually the case when pions are in- 
volved, then the extra pole does not come close to the real-energy axis, causing a broad structure 
in the scattering cross section, which often is not even a clear resonance. When threshold is closer 
to the confinement-spectrum ground state, then the extra pole comes closer to the real-energy 
axis, but the real part of the pole position is below threshold, causing structure at threshold. 
When threshold is sufficiently close to the confinement-spectrum ground state, then the extra 
poles end up on the real-energy axis below threshold, where they give rise to bound states, or 
possibly virtual bound states close to threshold. In Fig. [7| we show the various possibilities. 

In order to study the pole displacements as a function of the coupling, we normalize A such 
that A = 1 corresponds to the physical situation in Fig. The confinement spectrum is chosen 
to have a ground state at 1.3 GeV and a level spacing of 0.4 GeV. We study meson-meson 
scattering in S wave. One of the two mesons has a fixed mass of 0.5 GeV, the other is varied in 
order to obtain different values for threshold. A first observation from Fig. [7| is that the extra 
pole resides at minus infinity imaginary energy when the meson-meson system is uncoupled from 
the confinement states. 

From Fig. [7^, where the second meson mass is equal to 0.2 GeV, so threshold is at 0.7 GeV, 
we observe that the extra pole comes out at (0.711 — i0.20) GeV for A = 1. In the cross section 
this gives rise to a very broad resonance- like structure. 

In Fig. [7|d we choose 0.31 GeV for the second meson mass, so threshold lies at 0.81 GeV. In 
this case we find the physical pole on the real-energy axis, some 73 MeV below threshold, but 
still before having reached threshold. Such a pole represents a virtual bound state, which reflects 
itself in some structure at threshold in the scattering cross section. Only for A = 1.2 the pole 
arrives at threshold. Notice also that the pole reaches the real-energy axis at 0.70 GeV, which is 
110 MeV below threshold. Such a phenomenon only occurs for S- wave scattering. 

In Fig. [7^ the second meson mass equals 0.48 GeV, which brings threshold to 0.98 GeV. The 
situation near threshold is now more confusing, but we have enlarged that part of the figure in 
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Figure 7: The extra poles w.r.t. threshold. 

the inset. From the inset of Fig. [7^ we learn that the physical pole comes out some 10 MeV 
below threshold, but this time as a real bound state. Note that in this case the pole arrives at 
the real axis 230 MeV below threshold. 

In the next section we discuss well-known experimental facts supporting the existence of extra 
poles in the scattering amplitude for strong interactions. 

5 The scalar mesons 

As extra poles show u mainly in S"-wave scattering, it has to be expected that scalar mesonic 
resonances must be the domain to search for them [HJ. 

For Kit and tttt scattering, threshold is far below the ground states of the respective con- 
finement spectra us and {uu + dctj /y/2, which are at 1.389 GeV resp. 1.287 GeV in our model 
with the parameters of Pj. This is the situation comparable to Fig. [T^i. Consequently, the poles 
associated with the 1^(791) and the /o(400-600) resonances are still deep in the second Riemann 
sheet, causing broad structures in the scattering cross sections. In Ref. |B], the respective poles 
were found at (727— i 263) MeV and (470 — i 208) MeV, well explaining the scattering data. 

The KK threshold at some 0.99 GeV is much closer to the ss and (uu — dctj / y/2 ground 
states at 1.491 GeV and 1.287 GeV. The poles end up just below threshold, as bound states of 
the coupled systems [T4*lll7j . But ss also couples to tttc in our many-channel model jH], through 
the chain 
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ss — > KK — > — j= (uu + dctj — > tytt , (2) 

whereas (uu — dd) / \/2 also couples to in]. 

The ss channel couples to tttt through the chain (|2|L implying that the effective coupling 
is small fElGH]- Hence, the resulting pole displacement relative to the mentioned bound state 
below threshold, is small. In [Hj one finds (994— i 20) MeV for this pole, associated to the /o(980) 
resonance. 

In the case of the a (980) resonance, both KK and ttt] couple directly, through OZI-allowed 
pair creation/annihilation, to (uu — dd) j\[2. However, the coupling intensity for nrj is three 
times smaller than for KK. Consequently, we may deal with ttt] as a perturbative correction to 
KK here, resulting in an additional shift of the pole in the negative-imaginary direction in the 
complex-energy plane. In the full multi-channel calculation of Ref. jH], the pole corresponding to 
the a (980) was obtained at (968-2 28) MeV. 

6 The model's K matrix 

In its one but simplest form, the model's K matrix for low-energy elastic meson-meson scattering 
in I wave is given by the expression 



2a i X 2 jj l pjj(pa) J2 



\^mAa)\ 2 



Up) = -° m ~ Enl ; . (3) 

00 i (a) 

2a i X 2 fipj e (pa) n e (pa) J2 Wl T f 1 

E (P) - E nt c 

In formula (jSJ), p, fi, and E(p) respectively represent the linear momentum, reduced mass, and 
total invariant mass of the two-meson system. Furthermore, a stands for the average distance 
where quark-pair creation/annihilation takes place, J r n £ c is the radial part of the eigensolution 
of the confinement system with radial excitation n in i c wave and for eigenvalue E n £ c , j# and 
ri£ are the spherical Bessel resp. Neumann funtions, and A represents the coupling between the 
two-meson and the confinement systems. In Ref. ^2] we give a derivation of formula (j3J). 

One easily deduces for small A, assuming there exists a pole close to one of the confinement 
energy eigenvalues, say E v £ c , that the pole displacement, defined by 

A.^ c = E — E v i c , 

is given by the relation 
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-^»| 2 , (4) 
where 

pW i s 

the linear momentum corresponding to the energy E V £ C . 

Notice from relation that the imaginary part of the pole displacement is negative, as it 
should be when E vic is above threshold. 

It should also be clear at this stage that, by letting A increase, one may follow the pole's 
trajectory, but when next A is decreased, the pole nicely returns to E v g c . Hence, the extra poles 
do not show up this way. 

In the other limit of very large A, one ends up with a very simple relation for the pole positions, 
namely 

je(pa) + in e (pa) = , (5) 

which is the relation for infinite-hard-sphere scattering in i wave. 

For data analysis of meson-meson scattering data, the form of the K matrix in Eq. (jSJ) is as 
good as any Breit-Wigner expansion. But it has two advantages: 

• It automatically incorporates the extra poles, without any additional parameter. 

• In the limit A j 0, one obtains the confinement spectrum, equivalent to quenched-lattice 
spectra. 

In real data analysis, one may just take a few terms of the summation over the radial quantum 
number n, and approximate the rest of the summation by a constant. Since, furthermore, con- 
finement is not understood, the moduli squared of the eigenstates \ J r n £ c ( y a)\ 2 and the eigenvalues 
E n £ c turn into fit parameters for data analysis, to be adjusted to experiment. The resulting ex- 
perimentally determined confinement spectrum must coincide with the quenched-lattice spectra. 
This moreover implies that the light scalars are no issue for quenched-lattice calculations. Thus, 
a perfect mediator between experiment and low-energy QCD is created through expression ©• 

Finally, note that the expression (jSJ) could as well serve for baryon-meson [20] or baryon- 
antibaryon scattering. 

7 Conclusions 

Through unitarization, transition matrices can be constructed describing meson-meson scatter- 
ing. Here we discussed a method which reproduces low-energy data for those two-meson systems 
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that couple through OZI-allowed processes to nonexotic confined states. By studying its pole 
structure, we discovered that there exist two types of singularities in the scattering amplitude 
when analytically continued to complex energies. One type of poles can be directly related to the 
confinement spectrum, even for cases where pole displacements are large. In the limit of large 
coupling, these poles are equivalent to scattering from a hard sphere. Poles of the other type 
stem from the background, and are mainly important for scattering in S wave, the light scalar 
mesons being their most important manifestation, as well as the recently observed BABAR 
resonance [TUJ[T2]. 

When applied in data analysis, the model's K matrix might build the bridge between exper- 
iment and QCD lattice calculations. 
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